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g;; ■ Abstract 

^ \ The spin Hamiltonian that describes the molecular magnet Feg has biax- 

ial symmetry with mutually perpendicular easy, medium, and hard magnetic 
axes. Previous calculations of the ground state tunnel splittings in the pres- 
ence of a magnetic field along the hard axis are extended, and the meaning 

of the previously discovered oscillation of this splitting is further clarified. 
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I. INTRODUCTION 

The molecular magnetic cluster [(tacn)6Fe802(OH)i2]^"^ (henceforth abbreviated to Feg) 
is approximately described by a spin Hamiltonian of the type 



n = hJl + k^Jl - gfisJ ■ H, (1.1) 

where J is a dimensionless spin operator, ki > k2 > 0, and H is an externally applied 
magnetic field. For Feg, J = 10, ki ^ 0.33 K, and k2 ~ 0.22 K. The zero-field Hamiltonian 
has biaxial symmetry with easy, medium, and hard axes along x, y, and z respectively 0]. 

The magnetization dynamics of crystals of Feg and of a similar compound, Mni2-acetate, 
show some remarkable phenomena [0,^,0, which are not fully understood, and in this au- 
thor's view require a consideration of the spin-phonon interaction in addition to the pure 
spin Hamiltonian (|1.1| ) 0,|^. It is not the purpose of this paper to discuss these issues. 
It seems that at least some aspects of experimental behaviour can be understood on 
the basis of the spin Hamiltonian alone, and it is therefore worthwhile to address the sim- 
pler one-particle-quantum-mechanics problem of understanding the eigenvalue spectrum of 
Eq. ( |1 . 1|) . Several authors have taken this position, and that is our motivation in this paper 
too. 

In particular, we wish to focus on the behaviour of the two lowest energy levels when the 
external magnetic field H is along the hard axis, z. It was predicted in Ref. [^] that the 
tunnel splitting between these two levels would oscillate as a function of the field strength, 
and these oscillations appear to have been seen by Wernsdorfer and Sessoli IQQ. Very 
briefly, for H < He = 2kiJ/g^B, the spin system has two classical ground states, in which 
the spin is canted away from the ±x directions toward z. The classical degeneracy of these 
states is lifted by quantum tunneling. The tunnel splitting. A, oscillates as a function of H, 
vanishing at 2J field values lying in the interval {—H*, H*), where 

H* = {1- Xy/^Hc, (1.2) 

and A = /c2/fci- (This quantity, along with other important parameter combinations, is 
tabulated in Table I.) The vanishing fields are located symmetrically around if = 0; thus, 
in the case of half-integral J, if = is a vanishing point, in accord with Kramers' theorem. 
The above calculation was done using instanton methods, wherein the oscillation in A 
arises from the presence of two instantons, both of which have a complex Euclidean action. 
The imaginary parts of this action are equal and opposite for the two instantons (in a suitable 
gauge), and the vanishing of A can be understood as due to superposition or destructive 
interference between the corresponding amplitudes ||lT,12 . 



Recently, Chudnovsky and Hidalgo |]T3[ have examined this problem for H > H*, and 
noted that in this field range the splitting grows monotonically instead of oscillating. This 
is said by them to be a new topological effect or quantum interference phenomenon, arising 
from a new type of instanton. Indeed, their calculation of the splitting is quite intricate, and 
the imaginary part of the action vanishes due to a cancellation of two seemingly unrelated 
integrals, which they justifiably refer to as remarkable. 

In fact, the monotonicity of A for H > H* is already known and, we shall argue, entirely 
expected on physical grounds. This is done in Sec. II. We shall also present another way of 



calculating A, which is a simple extension of that in Ref. [|T0|, requiring nothing more than 
careful analytic continuation and contour integration. The imaginary part of the action 
vanishes in a simple and unremarkable way in this calculation. This is done in Sec. III. 

II. PHYSICAL REASON FOR OSCILLATION OF TUNNEL SPLITTING 



As noted in Ref. ||T^, there is another way to understand the vanishing of the tunnel 



splitting. For H || z, the Hamiltonian (|1 . 1|) becomes 

n = hJl + k^Jl - QfiBHj,, (2.1) 

which is invariant under a 180° rotation about z. This symmetry is reflected in the selection 
rule {m\T-C\m ±1) = 0, where \m) is a Jz eigenstate. The Hamiltonian divides into two 
disjoint subspaces, V+, spanned by m = J, J — 2, . . ., and V_, spanned by m = J — I, 
J — 3, .... The lowest eigenvalues (among others) from different spaces can cross as H is 
varied. The vanishing of the tunnel splitting is no more than such a level crossing. It was 
shown in Ref. [|14[ (see also Fig. 1 there) that there are exactly 2 J level crossings as H varies 
from —He to He- Since our previous instanton calculation ||T^ flnds just as many crossings 
in the smaller fleld range {—H*,H*), it follows that there cannot be any more crossings or 
oscillations in A for H > H* (or H < —H*). Chudnovsky and Hidalgo's findings merely 
refiect this already known fact. 

That A grows monotonically for large H is also to be expected on physical grounds. As 
the field grows and approaches He, the classical ground state orientations approach closer 
to z. The angle between them decreases, as does the energy barrier bewteen them. The 
default behavior in such a case is for the splitting to grow. It is the oscillation in A which 
is surprising, not the monotonic growth. 

The above argument shows that the vanishing of A is really not a topological effect in 
that it is not robust against perturbations such as a small misalignment of the magnetic field. 
Nevertheless, the instanton method provides an effective way of approximately calculating 
the splitting and the crossing points. 

III. INSTANTON CALCULATION OF SPLITTING FOR H > H* 

We now show how to adapt the calculation of the tunnel splitting A in Ref. |]TU1 , where we 



only gave explicit results for \H\ < H*, to the case H > H*. We refer readers to our previous 
paper for the calculational set-up. We also follow our previous notation, and summarize the 
important parameters in Table I. 

In the instanton method, the splitting is given by 

A = \£io,e-% (3.1) 

i 

where Ui is a quantity with dimensions of frequency or energy, and Si is the Euclidean 
action along a semiclassical path, or instanton, connecting the two states between which one 
is tunneling. The sum runs over all possible instantons. In simple problems there is only 



one instanton, but the converse is possible, especially in cases with high symmetry. This is 
the case for our problem, where there are two instantons. 

The dominant behavior of A is controlled by the action S, so we will focus on finding only 
this quantity. For our problem, symmetry guarantees that the prefactors Ui are identical 
for the two instantons involved, and the vanishing or nonvanishing of A hinges entirely on 
whether or not the actions Si are real or complex, and if the latter, what the relative phase 
of the corresponding amplitudes is. Thus there is no need to find the prefactors explicitly. 

The spin paths in question can be given in terms of spherical polar coordinates 
{0{r), 4>{t)) where r is an imaginary time. For a given path, the action is given by 

S = f (tJ{l - cos 6)^ + E{e, 0)) dr, (3.2) 

where the dot denotes a r-derivative, and E{6, 0) is the classical energy or expectation value 
of 7i in the spin coherent state that is maximally aligned along the direction {6,4>). For 
Eq. ( p.l|) , we get, up to an additive constant, 

E{e, 0) = kiJ^ [(cos 9 - cos ^o)^ + A sin^ 6 sin^ 0] , (3.3) 

where cos^o = ^o = H/Hc, and A = /c2/^i- Note that uq and A are both less than unity. 
This energy has minima at (^,0) = (6'o, 0) and {60,1^), and the additive constant has been 
adjusted to make £" = at these minima. 

The instanton paths are those for which the action is stationary, and which thus obey 
the Euler-Lagrange equations, 

dE 
iJ9 sin 9 = — — = 2XkiJ'^ sin^ 9 sin cos 0, (3-4) 

a0 

dE r -1 

iJ(psm9 = — —- = 2kiJ'^sm9 (cos^ — cos^o) — Acos6'sin^ . (3-5) 

o9 L J 

The boundary conditions are that the path approach the two minima as r ^ ±00. 

Along an instanton, E is conserved, so that the orbit without regard to its time depen- 
dence can be found purely by using energy conservation. This point is of great utility in 
calculating the instanton action, for if we arrange for E to equal zero, the action is given 



entirely by the first term in Eq. (3^), which can be written as an integral over 



5'instanton = ^^ / (^ ^ COs6'(0)) rf0. (3.6) 

In this integral, it is no longer necessary to regard as running over the original instanton 
contour given by 0(r). Any deformation of this contour is allowed as long as it does not 



encounter any singularities. This was the procedure followed in Ref. ||I0[. Setting ii^ = in 
Eq. ( |3.3D we obtain 

Mo + zAV^sin0F(0) 

cos^ = . 2 , , (3-7) 

1 — A sm 

where 

F(0) = (l-M^-AsinV)^/^ (3 8) 



It is important to note that which branch of the square root is taken depends on whether 
the instanton starts (at r = — cx)) at = or = vr. Let us consider the former case. It 
is obvious from symmetry that there are two instanton paths, {Oj-{t),(J)±{t)), which wind 
around z in opposite senses. Examination of Eqs. (|3.4|) and (|3.5| ) shows that we must take 
F((f)) > for both + and — instantons. The calculation is direct in the case mq < (1 — ^Y , 
which is equivalent to H < H*, for then, the integration contour for (p can be taken to run 
from to ±7r while keeping F{(j)) real. The integration in Eq. (|3.(j|) yields 



^4 



tJr Zt liJ'S 



where 



tJr- 




7(l-ng)(l-A)+nov/A' 



In 



S^ = J7r{l-uo/Vl-X). 



(3.9) 

(3.10) 
(3.11) 



The splitting is thus proportional to exp(— 5*^) cos(5'j), which oscillates as H is increased as 
discussed above. Note that Si = when uq = (1 — A)^/^. 

In the case when uq > (1 — A)^'^, i.e., when H > H*, one may guess the answer 
by reasoning that it must be analytic in the magnetic field, and that there should be no 
imaginary part to the action, since A grows with H as argued in Sec. II. And in fact, we 
notice that Si is exactly what one would obtain if the denominators of the arguments of 
the logarithms in Sr were multiplied by —1, and the branches of the logarithm interpreted 
suitably. This leads us to write 
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Uq 



Va-v/(i-«o)(i-a), 



(3.12) 



This expression is real when H > H*. In the other case, it can be made to agree with 
Eqs. (|3.1CI|) and (|3.11|) if we stipulate that for x < 0, Inx is to be taken as In \x\ =F iT^ for S±. 
It turns out that Eq. (|3.12|) is correct. 

Some readers may find the above argument too glib and reasoned after the fact. Is it 
possible to obtain it by honest calculation? The immediate problem is that when uq > 
(1 — A)^/^, F((f)) is imaginary for in the real interval A^ = (0i,vr — 0i), where sin^ 0i = 
(1 — uqY/X. One can proceed just as in the low field case by using Eqs. (|3.6|) and ( p.TD , 
keeping real, cos 6* is then purely real in the interval Ai, and the resulting imaginary 
contribution to the action cancels that from outside Ai. This cancellation has the same 
unsatisfactory and accidental character as that in Chudnovsky and Hidalgo's calculation 
p!3| . Further, one cannot decide which branch of F[(j)) should be used in Ai solely from 



energy conservation. 

Another procedure, which also entails simpler integrals, is to find the r-dependence of 
the instanton, and evaluate 5* as an integral over r via Eq. ( p.2|) . We record the relevant 
formulas for the case H < H* first. Only one instanton, say the + one, need be considered 
explicitly. Substituting Eq. (|3.7|) in Eq. (|3.5|) we obtain 



= u;osin0F(0), 



(3.13) 



where Uq = 2fci JA^/^. This equation can be integrated easily. Defining z = {1 — UqY^'^ujqt, 
and 

G{z) = {l-ul-XtSinh^zy/^, (3.14) 

we obtain 

cos^ = -{l-ul-\Y/^tanhz/G{z), or, (3.15) 

tanhz = -(l-n^)^/2cos0/F(0). (3.16) 

It is easily verified that 0^0, vr, as r ^ ±cxd. 
We also obtain 

^ = (1- «2)i/2(i _ ^2 _ xy/^ech z/G\z), (3.17) 

1-ul- X-uoXsech'^z-i{l+uo)X^^^il-ul- Xy^'^sechz 

l-cose = (1-uo) ^— 2 . (3.18) 

(1 — Mo)(l ~ -^j ^ ""oAtanh z 



When mq > (1 — A)^/^, it is clear that Eqs. (|3.13|-prTBD must continue to hold as general 



analytic relations since they arise from integration of the equations of motion which are 
unchanged. The only subtlety is in the assignment of branches to the various multivalued 
quantities and functions involved. Properly speaking, the integral in Eq. ( p.6| ) should be 
taken over a contour which is suitably indented around the branch points. It is clearly 
permissible to parametrize this contour in terms of any other variable, and one choice which 
is often convenient is to use the original time-like variable r, but allow it to become complex. 
With this point in mind we first make F{(f)) analytic by using a two-sheeted plane with a 
branch cut joining the points and vr — 0i. We then integrate Eq. (|3.13|) as 



(I -7,2U/2 

^ , ; ,, d^, (3.19) 

sm (pi' [(p) 

where the contour F is shown in Fig. 1. We leave it to the reader to verify that a consistent 
phase assignment can be made to F{(f)) in such a way that along this contour, we have 

argF(0) = f,, t ^ '^i' '^ ^ "r '^'' (3-20) 

^ ^^^ l7r/2 0i<0<7r-0i. ^ ^ 

The corresponding contour C traced out in the z plane is drawn in Fig. 2. The easiest way to 
see this is that for outside the cut, the integrand in Eq. ( |3.19| ) is real, so the corresponding 
parts of C must be parallel to the real axis. Secondly, the change in 2; as runs across the 
cut is given by 

A. = e-/'p . ^,.f:^\.^,.,. ^'P=-'-- (3.21) 

J(pi sm0(sm — sm 0ij^/^ 

Using the arbitrary constant of integration to adjust the overall vertical position of C, we 
arrive at Fig. 2. We do not bother to ask how C is indented around the poles of tanh^; at 

6 



±27r/2, because it turns out that the integrand for the action integral is analytic at these 
points. 

The problem is now reduced to integrating the product of the factors (1 — cos 6') and 
dcjy/dz over the contour C. The factors are given by Eqs. ( p.l7|) and (p.l8|) , with the proviso 
that 

{l-ul-Xf/^ ^i{ul + \-lf'^. (3.22) 

The correctness of this choice follows by noting that for (j) = tc/2, (f) = uJo(l — Uq — \Y'^, 
which must be positive and imaginary given our contour C. Collecting together Eqs. ( p.22| ), 
( KT% , ( Fl8l) , and (U), we obtain 



I+Uq 



Ua + X-1 +uoXsech'^z- X^/'^{ul + X- 1)^/^(1 + Mo)sech 2 
c ((1 -ul){l - A) -ugAtanh^^) (1 - ug - Atanh^ z) 



X / "° ^/ ~ ^ ^ "°""""^ "~" K-o^^-^J' y.^ao...... ^^^^^ ^^_ ^323) 



All square roots in this formula are regarded as positive. We now introduce the definitions 

,9 (1 — Un)(l — A) ,9 1 — Mn 

tanh^zi = - ^^ -, tanh2z2 = — r^, 3.24 

UqX a 

and the attendant results 

sech Z2 = uosech zi = X'^^^iul + A - 1)^/^ (3.25) 

The result for 5*+ simplifies greatly in terms of these quantities. The apparent singularities 
at z = ±2:1 and z = ±Z2 are cancelled by factors in the numerator, and we end up with 

S^ ^ j {^-ul)nul + X-lf/^ I ^chz ^^_ 

MoA(l + Uq) Jc (sech z + sech zi)(sech z + sech Z2) 

The integrand now has no singularities in the strip — vr < Im z < n. The contour C may 
thus be deformed into the real line, which proves that S+ is purely real. The actual integral 
itself is elementary, and may be done with the help of Ref . |T5[ . The final result is given by 
Eq. ( p.l2|) . The same result is obtained for S_. 
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FIGURES 

FIG. 1. The integration contour F in the plane. The heavy hne denotes a branch cut, and 
the dashed section of the contour hes on the second Riemann sheet. 

FIG. 2. The contour C in the complex z plane. 



TABLES 
TABLE I. Summary of important parameter combinations 



Quantity Formula 

He 2kiJ/gfiB 

Uq, COS^o H /He 

A k2/ki 

H* (1 - xyi^Hc 

ujo 2j{kik2y/^ 

sin 01 [(l-n2)/A]V2 

z (1 - ng)i/2wor 
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Fig. 2 



